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ABSTRACT 

Both global dynamics and turbulence in magnetized weakly collisional cosmic plasmas 
are described by general magnetofluid equations that contain pressure anisotropics and 
heat fluxes that must be calculated from microscopic plasma kinetic theory. It is shown 
that even without a detailed calculation of the pressure anisotropy or the heat fluxes, 
one finds the macroscale dynamics to be gcnerically unstable to microscale Alfvenically 
polarized fluctuations. Two instabilities that can be treated this way are considered in 
detail: the parallel firehose instability (including the finite-Larmor-radius effects that 
determine the growth rate and scale of the fastest growing mode) and the gyrothermal 
instability (GTI). The latter is a new result — it is shown that a parallel ion heat 
fiux destabilizes Alfvenically polarized fluctuations even in the absence of the nega- 
tive pressure anisotropy required for the flrehose. The main physical conclusion is that 
both pressure anisotropics and heat fluxes associated with the macroscale dynamics 
trigger plasma microinstabilities and, therefore, their values will likely be set by the 
nonlinear evolution of these instabilities. Ideas for understanding this nonlinear evo- 
lution are discussed. It is argued that cosmic plasmas will generically be "three-scale 
systems," comprising global dynamics, mesoscale turbulence and microscale plasma 
fluctuations. The astrophysical example of cool cores of galaxy clusters is considered 
quantitatively and it is noted that observations point to turbulence in clusters (veloc- 
ity, magnetic and temperature fluctuations) being in a marginal state with respect to 
plasma microinstabilities and so it is the plasma microphysics that is likely to set the 
heating and conduction properties of the intracluster medium. In particular, a lower 
bound on the scale of temperature fluctuations implied by the GTI is derived. 

Key v^rords: instabilities — magnetic fields — MHD — plasmas — turbulence — galaxies: 
clusters: general. 



1 INTRODUCTION 

Many astrophysical plasmas are not sufficiently colli- 
sional to be described by the standard fl uid equations 
of m aKnetohydrodynami cs (MHD) (see, e.g., 
ISche kochihin & Cowle j l2006l : ISharma et all 



Balbud 12004 
20061 . boOTi ). 



When the collision frequency u is smaller than the Lar- 
mor frequency n = eB/mc of the particle gyration about 
the magnetic-field lines, the plasma becomes magnetized: 
pressure and heat flux are now tensors that depend on 
the local direction of the magnetic field. This complica- 
tion leads to three significant physical effects. Firstly, on 
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the macroscopic scales, the momentum and heat transport 
become highly anisotropic with respect to the magnetic- 
field direction. Secondly, old MHD instabilities, like the 
MRI, th at are believe d to ex cite turbulence in astrophysical 
systems (iBalbus fc Hawlcv 1 9981). are significantly modified 
dQuataert et all 12002. : .Sharma et all [20031; llsla m fc BalbusI 
l200 5|) and n ew one s appear: MT I Ifialbus 2000l). MVI 
(Bai bu3[200i ). HBI llQuataertlbOOl ). Thirdly, a host of su- 
perfast microscale instabilities exi st that are directly driven 
by the pressure an isotropies (see ISchekochihin etaL I l2005l : 
[Sh arma et al.ll2006l . and references therein) and, as we are 
about to discover, also by heat fluxes. 

The presence of microscale instabilities especially opens 
a fundamental problem: the equations one tends to use to 
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describe the macroscopic dynamics of magnetized plasma, 
be they fluid or kinetic, are derived in the long-wavelength 
limit (fcp <C 1, where p is the Larmor radius; see iKulsrudI 
1 19831 ) and turn out to be ill-posed because in this limit 
the microinstabilities hav e growth rates proportional to k 
jSchekochihin et al.|[2005l l. In order to regularize them at 
small scales, one has to take into account eff'ects associated 
with the finite Larmor radius (FLR), which requires fairly 
complicated kinetic theory and typically means that the full 
multiscale problem is analytically hard and numerically in- 
tractable. Ideally, one would like to have an effective mean- 
field theory, with the microscale fiuctuations analytically av- 
eraged to produce some form of closure for the momentum 
and heat transport. This has not been achieved yet, but an 
educated guess about the form of such a closure can me 
made, based on the idea that the system should always find 
itself in the marginal state with respect to the microinsta- 
bilities JSharma et all |2006|. l2007l: ISchekochihin fc CowlevI 



l2006l : lLvutikovll2007l : lKunz et aLlbOlol ) 



In this paper, we attempt to make progress in setting up 
the theoretical framework for astrophysical plasma dynam- 
ics by addressing three basic questions: what is the general 
form of the dynamical equations that we are attempting 
to approximate? what can be learned about the microin- 
stabilities under the most general assumptions? what con- 
straints do their marginal stability conditions impose on the 
allowed macroscopic states of the plasma? The first of these 
questions is addressed in section [2] the second in section [3l 
where an old (firehose) and a new (gyrothermal) instabilities 
of Alfvenically polarized perturbations are derived. Possible 
ways of thinking about the nonlinear physics of these mi- 
croinstabilities are proposed in section [4] The physical con- 
clusions are summarized in section O including a discussion 
of the relevance of all this in galaxy cluster cores (as a case 
study of a multiscale astrophysical plasma system). 



2 EQUATIONS FOR PLASMA DYNAMICS 

Let us consider a two-species fully ionized plasma. In the 
completely general case (assuming only quasineutrality) , the 
evolution of ion density n and fiow velocity u is governed by 
the following equations 

An 
dt 



= — nV ■ u, 
Au 



mn —— = — ^ 



At 
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where m is the ion mass. A/ At — d/dt + u- V the convective 
derivative, I the unit dyadic, B the magnetic field and P 
the plasma pressure tensor. It is via P that all the kinetic 
physics comes in: in general, P is the sum of the ion and elec- 
tron pressures and for each species, it is P = m J d^vvvf, 
calculated from the distribution function f{t,r,v), which is 
the solution of the kinetic equation for that species. Note 
that V is the peculiar velocity, i.e., the particle's velocity in 
a frame moving with the mean fiow velocity u. 

Thus, the challenge is to calculate P. This typically in- 
volves setting up an asymptotic expansion of the kinetic 
equation with respect to one or several of the small param- 
eters available for the plasma under the set of macroscopic 
conditions of interest. Many such expansions for magnetized 



plasma exi st, correspondin g to various physical regimes : 
collis i onal (iBraginskiil 19651: jMikhailovsku fc TsvpinI Il971 , 



1 1984 ICatto fc Simakoy 2004), long- wav e length collision - 
less, or drift-kinetic (|Chew et al.1 1 19561 : iKulsrudI Il983l) 



short-wavelength anisotropic, or gyrokinetic ( Howes et al.l 
l2006l : ISchekochihin et al.ll2()()9l . and references therein), and 
more specialized versions of the above, appropriate for the 
treat ment of pressure-a nisotropy-drivcn instabilities: fire- 
hose JSchskochihin ,2008: Rosin et al., ,2010.) and mir- 



ror (Califano et al.l 200 j : Istomin et al.l 20091 : iRincon et al.l 
l2O10i r We do not at the moment wish to pick any one of 
these, but simply notice that in all of them, the equilibrium 
distribution function invariably turns out to be gyrotropic, 
i.e., independent of the phase angle of the particle's Lar- 
mor gyration. The only assumptions needed for that is that 
the characteristic frequencies uj for the evolution both of 
the equilibrium and of the perturbations thereof should be 
smaller than the ion Larmor frequency Q and the length 
scales of the equilibrium longer than the ion Larmor ra- 
dius p. If the pressure tensor is assumed to be determined 
purely by the gyrotropic lowest-order distribution, then it 
reduces to a diagonal form, P = p_l(I — bb) + p^^bb, where 
b — B/B, and the perpendicular and parallel pressures are 
p± — mj A^v {v^/2)f and = m j d^uuy/. These pres- 
sures can be shown to satisfy th e so-called CGL equations: 
for each p article species, they are (|Chew et al.ll 19561 : Kulsrudl 
Il983,: .Snyder et al.lll997l : see Appendix[B]for a simple deriva- 
tion) 



At nB 
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where u is the collision frequency, q± = m J A'^vv\\{v^/2)f 
and = m J A'^v v^^f are the parallel fiuxes of the perpen- 
dicular and parallel heat and Qj^ — bq±, qy — bgy. 

As mentioned above, pressure anisotropics p± ^ P\\ 7^ 
lead to instabilities whose peak growth rates occur at scales 
smaller than those allowed by the validity of the diagonal ap- 
proxi mation for P and are not captured by this approxima- 
tion jSchekochihin et al.ll2005l ). The instabilities are regular- 
ized by the FLR effects, so it is natural to resort to FLR cor- 
rections in the plasma pressure tensor ([Snvder fc HammettI 
I2OOII : iRamosllioOsI : iPassot fc Sulemll2007l ). To lowest order 
in Lj/il and kpi, this is quite easy to do and the result, a 
simple derivation of which is given in Appendix [X] is 

P = pxl-(px-P||)66 + G, (5) 
G = ^ [b X S • (I 36b) - (I 366) -5 x 6] 

+ i[6(ax6) + (ax6)6], (6) 
where the auxiliary tensor S and vector a are 

S = (px Vu + Vq^) + (pxVu + Vq^f , (7) 

a- = ip^-pii) (^^ + 6- Vn^ +(3gx-g[|)6- V6. (8) 

Each plasma species contributes a pressure tensor of the 
form ((5]). In general, electron pressures are comparable to 
ion pressures, but it is not hard to show that the electrons' 
contribution to the FLR term G is smaller than the ions' by 

a factor of {rrie/mi)^''^ . 
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Note that if one sets p± — P\\ =0 and 3q± — =0 
(as would be the case for an isotropic equihbrium distri- 
bution and colUsional heat fluxes), the FLR term G in 
equation ((5]) is readily recognized as the so called "gyro- 
viscosity^ tensor, first obtained (in the collisional limit) by 
iBraginskiH l|l965h (he assumed sonic flows and found just 
the Vtx terms; the heat flux terms w ere introduced later by 
iMikhailovskii fc TsvpinI (|l97ll . Il984l ) to accommodate sub- 
sonic flows). 

Thus, the momentum equation ([2]) has the form 



du 

mn —— = — V 



dt 



+ 8^) 



+ V 



(9) 



where G is given by equation ((6)1 . Now we need an evolution 
equation for the magnetic field. Faraday's law reads 



dB 

dt 



-cV X E, 



(10) 



where E is the electric field. The electron momentum equa- 
tion is used to calculated E. Since the electron mass is small 
compared to the ion mass, to lowest order in {me/mi)^^^ this 
reduces to the force balance 



E + 



Ue X B 



0. 



(11) 



The electron density Ue is related to the ion density n by the 
quasineutrality of the plasma. Tie = Zn (the ion charge is 
Z times electron charge e). The electron flow velocity Ue is 
related to the ion flow velocity ti by ite = it — j'/erie, where, 
using Ampere's law, the current density is j — cV x B/4-k. 
Finally, since the FLR terms in the electron pressure tensor 
are negligible to lowest order in (me/mi)^/^, we have = 
pxel — (p±e — P||e)bb. Assembling all this together, we gelQ 
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Note that c/eue = B/mnQ., where m, n and Q. axe ion mass, 
density and Larmor frequency, respectively. 

We will not be preoccupied here with the determination 
of the pressures and heat fluxes (which is necessary to close 
the set of equations we have written down). Depending on 
the physical regime one is interest ed in, they can either be 
calculated in the collisional limit l|Braginskiil 1 19651 ) or Lan- 
dau fluid closures ca n be devised for them, appropriate for a 
collis i onless plasma (ISnvder et al.ll 19971: Snyder fc HammettI 
I2OOII : iRamodlioOsI : iPassot fc Sulemll2007l ). Instead of wad- 
ing into this rather complex subject, we will inquire what 
can be learned just from the general form of the equations 
of plasma dynamics outlined above. 



^ See Appendix [C] for the demonstration that this equation con- 
serves magnetic flux except at very smaU scales, where electron 
pressure anisotropy can lead to violation of flux freezing. 



3 FIREHOSE AND GYROTHERMAL 
INSTABILITIES 

In any given astrophysical problem, one might find some 
macroscale solution of the equations of section [S] describ- 
ing the large-scale dynamics. Such solutions turn out to be 
generically unstable to perturbations with large wavenum- 
bers and high frequencies (much larger than the fiuid 
turnover rates uj ^ | Vu|). In general, showing this involves 
having to perturb all quantities, including the pressures and 
the heat fluxes, which requires a kinetic closure. However, 
there is a class of perturbations whose stability does not 
depend on the details of kinetic theory. 

Let us start by perturbing the momentum equation (|9]) . 
We assume the perturbation to be oc exp(— itjt-|-ifc-r). In our 
perturbation theory, we will always consider terms contain- 
ing u) and fc to be dominant in comparison with the terms 
containing time derivatives or gradients of the macroscale 
quantities. Thus, from equation we get, noting that 
V ■ B = implies V ■ b = ■ VB/B, 



mnujSu — k± | Sp± + 

'fell I pi_ -pii + 



B5B\ ^, 
I +bk 
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Note that 5S = ip±{k5V + 5Vk), where 5V - 
b5q±)/p±. Therefore, from equation (|6}, 



A:||b+ (fex X b) - (fcx X b) ('fc||b+ ^ 



Su+ {q±Sb + 



+ — [fcii (Ser X b) + fcx ■ {Scr x b) b] 



■5V 



(15) 



Scr = — i [{p± — p\\) [ujSb — — (3gx — q\\)k\\?>\)\ . (16) 

In the above equations, SB = 5B<^\ and <5b = SB ^/ B, where 
SB satisfies the perturbed equation (|12|l : 

u = — fc||5ux + b (fcx ■ <5ux) 
B 



+ 



ik« 



mn Q 
+ (fcx X b) 



+ ^ ) (fe >^3b) 



(17) 
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B^\ SB 



Examining equations (|14H17[) . we observe that in the 
simplest case of fcx = 0, the Alfvenically polarized pertur- 
bations decouple from the compressive/slow- wave-polarized 
perturbations {Sn, SB, Sp±, (Jpy, Sq± and 5g||). No ki- 
netic physics is required to study the stability of Alfvenic 
perturbations, which satisfy 



mnujSui 



ik« 



B \ 

-P\\i+P±c-P\\e+ ^] Sb (18) 



P\^^Su± + (pxi - P\\i) -rr^b" (2'7xi - 1\\i)Sb 



ujSb = — fcii 5u± -f 



ikl 



'Plle + ^]{bxSb), (19) 



mn Q, 

where we have restored species indices on pressures and heat 
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fluxes; note that only ion FLR terms are kept in equa- 
tion (|18p . In the absence of FLR efi'ects, equations (|18p 
and (|19p describe Alfven waves with propagation speed 
modifled by the pressure anisotropy. When — py < 
-5^47!-, it gives rise to the well known firehose insta- 
bility with a growth rate 7 oc fc y jRosenbluthl 1956 ; 
Chandrasekhar. Kaufman fc Watsor] Il958l : iParkerl 195g 



Vedenov fc 1958l l. The FLR gives rise to a dis- 



persive correcti on that sets the wa.venumber of the fastest- 
grow ing mode (| Kennel fc SagdeevI Il967l : iDavidson k, Vollj 
1 19681 ). but it also contains a contribution from the heat 
fluxes, which lead to a new instability. 

Let us combine equations (|18|) and (|19|) and non- 
dimensionalize everything: 

2 V P 



5b + 



[(1 - S)uj + fcrr] {b X Sb) , 



(20) 

where uj = ui/Q, k — k\\p, p = Uth/f^, '^th = (2p\\ilmn)^l'^ . 
The problem has four physical dimensionless parameters 
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(21) 
(22) 



but, in fact, only two matter because only the combination 
A + 2/j3 figures in equation (|20|) and 5 will turn out not to 
be of much consequence. The resulting dispersion relation is 



2 fc / A 2 
u A + - 

2 V /3 



^[(l-5)a; + fcrT 



(23) 



This has four roots of which two can be unstable: 
CD = ±^(1-5) 
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i\k\ 

7f 



A + - I =F fcFT 



fc2 



(24) 



(we will henceforth refer to the positive/negative frequency 
modes as "+/— modes"). The instability occurs for k such 
that the expression under the square root is positive. De- 
manding that the interval of such wavenumbers is non-empty 
gives the necessary and sufficient condition for instability: 



Ai 



A + 



> 0. 



(25) 



3.1 Firehose instability 

We observe first that if the heat fiuxes are negligible, ^ 
\A + 2/j3\, this condition is satisfied for A + 2/13 < and 
we have the standard parallel (k ± = 0) firehos e dispersion 
relation (Kcrmcl & SaEdcov 19671: [Davidson fc Volk 1968): 
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where fco is the cutoff wavenumber and each of the -|- and — 
modes has two peaks of the growth rate occurring symmet- 
rically at kp — ±ko/V2 (see figure [T]a), where 
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Figure 1. Growth rates of the firehose and gyrothermal insta- 
bilities (equation i24l ) in three qualitatively different regimes: 
(a) pure firehose, Ft = 0; (b) GTI combined with firehose, 
A + 2//3 < 0; (c) pure GTI, A -I- 2//3 > (firehose stable). We 
have set S = A — 2//3, i.e., p±f, — P\\s = 0. 



A-h 



(28) 



Note that here and everywhere else, we assume that A is 
not too close to 1. 



3.2 Gyrothermal instability 

The situation becomes more complicated when the heat 
fiuxes are not negligible. Let us assume, without loss of gen- 
erality, that Ft > (otherwise, change the sign of the par- 
allel spatial coordinate). There are two unstable intervals: 



4 ( Ft -f %/A 
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: (Ft - Va) 
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(29) 



(30) 



When A + 2//3 ^ 0, these intervals intersect and contain k — 
0, otherwise they are disjoint (see figure 0),c). Computing 
their peak growth rates and corresponding wavenumbers is 
straightforward. Here we consider two interesting limits. 

When Ft » |A-f2//3|, we have, for the + and — modes 
respectively: 



fcp = =F 



6F2 



(1-5)2 



Tmax 



513 ■ 



(31) 



We see that an instability is present that is driven purely by 
heat fiuxes, even when the pressure anisotropy is neutralized 
by the tension force (A = —2/13). This is the purest form 
of the gyrothermal instability (GTI), which, as far as we 
know, has not been previously reported in the literature. In 
the more general case when the pressure anisotropy is not 
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negligible, the GTI operates in conjunction with the firehose. 
The condition (|25p means that GTI can be operative even 
when A + 2//3 > 0, a regime in which the Alfven waves have 
previously been believed to be stable. 

The second important limit is the case when GTI is 
close to marginal stability, A — )■ +0 (we are assuming that 
is finite, so the firehose is stable in this limit). According 
to equations (|29l) and (|30|) . the instability intervals in this 
limit shrink to the immediate vicinity of just two wavenum- 
bers: 



A_ 

1 J, 



where the upper sign is for the + mode, the lower for the 
— mode. This is a very different behaviour from the fire- 
hose, for which the interval of growing modes moves to ever 
longer wavelengths as marginal stability (A + 2//3 — >■ —0) is 
approached (see equation (|27|) ). i.e., the firehose stops being 
a microscale instability in this limit. In contrast, the GTI al- 
ways excites Alfvenic fluctuations at very short wavelengths. 

Finally, we note that the assumption in our derivation 
that uj/Q. <Si 1 and fcyp ^ 1 imposes constraints on the val- 
ues of our dimensionless parameters that we are allowed to 
consider: |A -|-2//3| ^ 1 for the flrehose and Ft <Si 1 for the 
GTI. The expressions for maximum growth rates and corre- 
sponding wavenumbers derived above (equations (|28p . (|31|) 
and (|32l) ') provide guidance on the relative smallness of all 
these quantities and, therefore, on the ordering schemes that 
can be pursued in weakly nonlinear theorie s (one example 
is the ordering adopted bv lRosin et"alll201(]| ) . 



4Ft/A 



(32) 



4 NONLINEAR EVOLUTION 

Nonlinear theories of pressure-anisotropy-driven plasma 
instabilities are in their infancy, but most of them agree 
that the net result is to drive th e anisotropics towards 
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). Observational evidence from 



the solar wind strongly p 

jKasper. Lazarus &: Garvl I2OO2I: 



oints i n the same d ir ection 
Hellinger et"ai] I2OO6I : 



iMatteini et al.ll2007l : lBale et al.ll200"9f ) 

If we assume that this is what happens in the case of 
the firehose and gyrothermal instabilities, then the marginal 
state A = (see equation (|25|) ) implies a certain relation- 
ship between the heat fiuxes and the pressure anisotropy 
in the nonlinear regime. In order to find the way in which 
the system contrives to set up this relationship, we must 
first examine the physical mechanisms that determine A, 
q± and gy . 

Subtracting equation ^ from equation ((21, we get 



d(p± - P\\ 
&t 



(px + 2p|| 
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(3P|| 



Q± - qii ) - 3gx V ■ h - iv{pA 
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n dt 

-P\\) 
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This tells us that there are three sources of pressure 
anisotropy: changing magnetic-field strength (changes in p± 
have to match changes in B to maintain conservation of the 
first adiabatic invariant for each particle, /i = mv\/2B), 
compression/rarefaction, and heat fiuxes. 



If we assume for a moment that the collision rate is 
larger than the rate of change of all fields, then the differ- 
ences between and py in equation (|33|) can be neglected 
everywhere except the coUisional term and so the steady- 
state average pressure anisotropy satisfies 



A=i 1^ 
v\BAt 



2 1 An V ■ [b(gx - gy)] -^3gxV ■ b 

3 n df 3p|| 



(34) 

Note that if we use equations \12\ and ([l| (neglecting FLR 
terms in the induction equation) to express the rates of 
change of B and n in the right-hand side o f equation (|34p . 
the first two terms are the lBraginskiil (|l965') parallel viscous 
stress. T he last term is the heat-flux co rrection to it intro- 
duced bv lMikhailovskii fc Tsvpiiil l|l97ll . Il984l ) for subsonic 
flows. Under the same assumption of high coUisionality, the 
heat fluxes ar^ 

gx = ig||=-in^b.Vr, (35) 

where T = p/n and p = (2/3)px + (1/3)P|| • 

As we showed in section |3] the slow macroscale motions 
that produce this A and these heat fluxes are unstable to 
microscale perturbations, in particular, the Al f venic ones ex- 
cited by the flrehose/GTI. ISchekochihin et all l|200^ showed 
that the way a sea of small-scale Alfvenic fluctuations can 
change a large-scale driven anisotropy is by growing secu- 
larly with time and thus producing a finite change in the 
average field strength: 



1 AB 
B~At 



J_ABo 1 9|5b|2 
B^~dr ^ 2 dt '■ 



(36) 



where the overbar denotes a small-scale average. Bo is the 
slowly changing macroscale held and Sb = 5B±/Bo is the 
fast microscale Alfvenic perturbation of it. Let us replace 
the magnetic term in equation (|34[) with its average given 
by equation (|36p . Even though the fluctuation amplitude 
is small, the nonlinear feedback will produce a finite con- 
tribution to A if the fluctuation energy grows secularly, 
\Sb\^ ~ 7ot, where 70 is the typical rate of change of Bq. 
There does not appear to be any other way for the small 
Alfvenic fluctuations to affect the average macroscopic pres- 
sure anisotropy or heat fluxes. 

In the case of the pure firehose instability (no heat 
fluxes), the nonlinear feedback described above cancels the 
negative pressure anisotropy that triggered the flrehose and 
pushes the system towards A + 2/13 — >■ —0. If heat fluxes are 
present, the marginal state of the GTI requires A-I-2//3 > 



2 The numerical prefactor in the last expression in equation l|35|l 
depends on the exact form of the collision operator used and 
is not relevant to our discussion. The same applies to numer- 
ical coefficients in equation 1134 l l and so, to preserve consis- 
tency, we have given the va lues obtained by using the Lorentz 
operator llRosin et al.ll2010h . The more precise coefficients for 
10ns are 25/32 in equation 3075/1068 in front of the 

first two terms in equation 1341 1. and 1823/1068 in front of 
the heat flux terms in the same equation; the ion collision fre- 
quency is 1^ = 4\/7rne^A/3 m^/2T3''2, where A is the Coulomb 
logarithm (B raginskii 19651: iMikhailovskii fc Tsvpinlll97ll , Il984l: 
ICatto fc S'imakovll2004l) . 
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(A — )■ +0; see equation (|25|) ). This can still be achieved by 
secularly growing Alfvenic fluctuations (which, unlike for the 
firehose, now have a definite scale unaffec ted by the pres sure 
anisotropy; this is explored further in Rosin et al.lboiol '). 

A remarkable consequence of this predicted tendency 
for a system to develop positive pressure anisotropy to cancel 
the destabilizing effect of heat fluxes is that instabilities as- 
sociated with A > (e.g., mirror) could perhaps be triggered 
as secondary instabilities of the saturated state of the GTI. 
One might imagine a sea of Alfvenic fluctuations attempting 
to neutralize the GTI and exciting unstable mirror modes 
— this is feasible if the pressure anisotropy corresponding 
to the marginal state of the GTI exceeds the mirror sta- 
bility threshold: A ~ 2T'^ - 2//3 > l/P, i.e., F^, > 3/2/3. 
The mirror mode near its threshold A — — )■ +0 is 
polarized as a highly oblique slow wave: it has 5u\\ and 
with fciip ~ A - < k^p ~ (A - (see, 
e.g.. lHelling"eill2007l ). This suggests a three-scale system: a 
macroscale equilibrium, the microscale Alfvenic foam with 
fc|| p ~ Ft ~ 1/ yfP (see equation ((32])) driven by the GTI and 
producing an average pressure anisotropy, and a mesoscale 
near-threshold mirror turbulence driven by that anisotropy 
and, because of scale separation, probably otherwise discon- 
nected from the Alfvenic modes. Finding out how they all 
coexist and how the mirror saturates requires a systematic 
kinetic calculation, which will be attempted elsewhere. 

Finally, as an alternative to the above considerations, 
we should perhaps mention the possibility of strong non- 
linear distortions of the magnetic fleld (5b ~ 1) that could 
reorient the field so as to minimize the parallel ion temper- 
ature gradient and thus switch off or weaken the GTI — 
on large scales, such behaviour has been observed in sim- 
ulations of another, macroscale, instability driven by the 
parallel (electron) heat fiux and buoyancy force, called the 
hcat-fiux-buoyaiicy inst ability, or HBI (Shar ma et al.,.2009i : 
IParrish et ai.1i2009i : iBogdanovic et al.ll2009l ) . 



5 PHYSICAL AND ASTROPHYSICAL 
CONSIDERATIONS 

5.1 Physical conclusions 

The main physical conclusion is that parallel heat fluxes 
can directly drive microscale instabilities in magnetized as- 
trophysical plasmas. This can happen in two ways. 

First, as follows from equation (|33p . plasma pressure 
anisotropy can be driven by heat fluxes, so flrehose, mirror 
and the rest of the microinstabilities due to p±_ — P\\ 7^ 
can be triggered not just by plasma motions, but also by 
parallel temperature gradients. Although perhaps not much 
discussed explicitly, this instability mechanism is not partic- 
ularly surprising and it is implicitly present in the existing 
analytical and numeric al models based on CGL equations 
with heat fluxes (e.g., Snyder et al.l Il997l : iQuataert et al.l 
l2002l : ISharma et al.ll2006l . l2007l ;i. 

A more interesting and, we believe, novel instability 



importantly, the GTI persists even when the flrehose is sta- 
ble, so the firehose marginal stability condition has to be 
replaced by the GTI marginal stability condition involving 
both the pressure anisotropy and the ion heat fiux (equa- 
tion pS)) ). 

The GTI is distinct from the two other in- 
stabilities associated with the presence of tempera- 
ture gradients and rece ntly ex p lored in astrophysical 
contexts — the MTI ijBalbud l2000l : IParrish fc Stond 
20071: IParrish et al I l2008l) and the HB I jQuataerd [iooi: 



Sharma et al. ' '2009'; 'Pa rrish et all l2009l : IBogdanovic et all 
2009; Ruszkowski & OlJ l2009l l. The latter are driven by 
buoyancy and are essentially mac roscale fiui d instabilities , 
hke MRI jBalbus fc Hawlevlll998l ) or MVI l|Balbusll20oi '). 
They are also much slower than the GTI, which is a mi- 
croscale plasma instability belonging to the same class as 
the flrehose, with peak growth rate a fraction of the cy- 
clotron frequency. Since such an instability can be triggered 
by the presence of a heat flux, one might wonder whether 
in the same way that large-scale pressure anisotropy could 
be conjectured always to be determined by the marginal 
stabi li ty co n ditions of the microinstabilities JSharma et al.l 
l 2006l. I2OO7I: ISchekochihin fc CowlevI I2OO6I : iLvutikovl 120071 : 
iKunz et alll2010l ). the heat fluxes as well should be con- 
strained by the marginal stability conditions of the GTI and, 
perhaps, other such instabilities. We stress, however, that, 
whereas this might be a reasonable interim course of action, 
it by no means excuses us from the task of finding out how 
GTI and the rest of the instabilities behave and saturate on 
the microphysical level (see discussion in section |4|. 



5.2 An astrophysical example: galaxy clusters 

A detailed development of applications to concrete astro- 
phy sical systems falls outside the scope of this paper (see, 
e.g. , iKunz et al.ll201ol ) . However, it is, perhaps, illuminating 
to provide a few estimates of the role the GTI might play 
in cool cores of relaxed galaxy clusters, a good example of 
a real astrophysical plasma for which a sufficient amount of 
observational evidence exists to enable a quantitative dis- 
cussion of the multiscale dynamics. 



5.2.1 Three-scale dynamics 

The conditions in the cluster cores are believed to 
be controlled by a balance between the radiative cool- 
ing and a reheating due perhaps to electron heat 
conduction from the bulk of the cluster and per- 
haps also to the t urbulence exc i ted by the active 
galac t ic nuclei (e.g., BinnevI 20031: Dennis fc ChandranI 
I2OO5!; IPeterson fc Fabian '2006"; McNamara fc N ulsenll2007l : 
,Guo et al. 2008; Ruszkowski fc Oh 2009 , and references 
therein). The plasma in the cores has the electron density 



mechanism is the destabilization of the Alfvenic perturba- 
tions by the ion parallel heat fluxes via the FLR effects in the 
plasma pressure tensor — we call this the gyrothermal in- 
stability (GTI). When the firehose is unstable, the GTI can 
substantially modify (increase) its growth rate, but more 



Jie in the range 10~ to 10~ cm~ at the radial distance 
of r ~ 10 kpc from the centre and about a factor of 10 less 
at the edge of the core at r ~ 100 kpc. The ion density is 
the same for a hydrogen plasma. The electron temperature 
Te is measured reasonably precisely and is of the order of 
a few keV, rising b y about a factor of 2 or 3 from r ~ 10 
kpc to 100 kpc (e.g . .lDavid et al.ll200ll:IVikhlinin et al.ll20"05l: 
iFabian et ai1l2006l : iteccardi fc Molendill2008l : ISanders et al.l 
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l2010al lb[). The ion temperature is not measured, but the ion- 
electron temperature equihbration turns out to be quite fast 
compared to all other relevant dynamics, so Ti ~ Te can 
reasonably be assumed. The unsolved macroscale problem 
is why the temperature does not drop lower in the centre — 
simple estimates suggest that the system should be vulnera- 
ble to a collapse onto the centre precipitated by the radiative 
cooling on a characteristic time scale of about 1 Gyr. 

This is where turbulent heat conductior(f| and turbu- 
lent heating are invoked as mechanisms that prevent the 
cooling catastrophe. The outer scale L of turbulent motions 
is believed to be between a few and a few tens of kpc, 
with corresponding velocities [/ of a few hundred kms~^ 
l|Enfihn fc Vogtlbood : ISanders et a"i]|2010d ). The turbulent 
motions lead to fluctuations in the magnetic-field strength 
and so excite pressure anisotropics, given by equation (I34|) . 
Equation (|12[) tells us that the typical rate of change of 
the field is comparable to the typical rate of strain ~ 
{U/ L)Ke^^'^ , where Re ~ ULv/v^^ is the Reynolds number 
(the maximum rate of strain that can affect the magnetic- 
field strength is at the viscous scale set by the parallel viscos- 
ity; see[Schckochihin & Cowlcv 2006 for a detailed explana- 
tion). Thus, we estimate the pressure anisotropy as follows: 



A-i^ReV^ 
1/ L 



1 uV^ U 



0.007 



ly L J vth 



( r 

Vo.Olcm-3/ 



IkeV 



1/4 



100 kms"^ 



10 kpc 



-1/2 



(37) 



where f is the ion collision rate. In view of the instability 
condition ((25|, whether this anisotropy will trigger plasma 
microinstabilities is decided by comparing it with 



P 



= 0.005 



B 



0.01 cm-» 



IkeV 



(38) 



The two numbers are remarkably close (obviously, only or- 
ders of magnitude matter here, given all the uncertain- 
ties). Thus the intracluster plasma teeters at the brink of 
marginal stability. In the unstable state, at the reference 
values B = 1 ^iG and Ti = 1 keV, the firehose (or GTI) will 
have growth times and peak-growth scales 



Tmax 
I -1 _ 



- (An,) 

^-1/2 



- 2 • lO'^s ~ 6hr, 
Pi ~ 700,000 km ~ 20npc 



(39) 
(40) 



(see section [3|. These are microscopic scales compared both 
to global cluster dynamics and intracluster turbulence. The 
implication is that the plasma instabilities should saturate 
and presumably contrive to return the intracluster medium 
to marginal stability instantaneously fast via an observation- 
ally invisible sea of nanoparsec-scale magnetic fluctuations. 



^ Since tlie cooling rate is oc rigTe and the relaxation rate 
of temperature gradients based on Spitzer conductivity is oc 
rig^T^^'^ l lSpitzerl 1 19621) , they cannot balance in a stable way, 
so Spitzer conduction by itself is not sufficient to explain the ab- 
sence of the cooling catastrophe. In contrast, turbulent heating 
controlled by the plasma instabilities via pressure anisotropy (as 
explained below) turns out to be a thermally stable mechanism 
for regulating cooling flows jKunz et al.ll201oli . 



Thus, a cluster core is a "three-scale system": global 
equilibrium profiles (10^ kpc, 10° Gyr) and turbulence 
(10^ kpc, 10^ Myr) constitute the macroscale magnetofluid 
dynamics of the intracluster mediumQ subject to trans- 
port properties controlled by "nanoscales" (10^ npc, 10^ hr), 
where plasma microinstabilities are excited. Their nonlinear 
behaviour sets the pressure anisotropy and probably also the 
heat fluxes. The pressure anisotropy determines the effective 
viscosity of the plasma and, therefore the heating rate; the 
heat fluxes determine the effective thermal conductivity — 
thus, neither the turbulence nor the global dynamics (e.g., 
temperature proflles for the cooling-core problem) can be 
computed correctly without a good theory or, at least, a 
good model prescription, for the effect of the microinstabili- 
ties on the macroscale dynamics. A similar three-scale situa- 
tion arises in most other weakly collisionalf] cosmic plasmas: 
e.g., accretion flows, solar wind, etc. 



5.2.2 Temperature fluctuations 

As we have shown in this paper, ion temperature gradi- 
ents, including ones due to temperature fluctuations, if they 
are there and if the associated parallel heat fluxes are large 
enough, will excite microinstabilities. The estimates of 7max 
and kp in section [5. 2. II still hold, by order of magnitude, for 
the GTI, so the instability is extremely fast and one should 
expect to find plasma close to the marginal state. We may es- 
timate (crudely), the minimum parallel temperature length 
scale allowed by the instability condition H25|) by requiring 
< 2//3 for stability and using equation (|35|) for the heat 
fluxes: 



It > 0.3 



0.01 cm- 



-1/2 



IkeV 



5/2 



B 



kpc, (41) 



where l^^ — b- V In T is the temperature scale. Note the very 
strong temperature dependence of this lower bound: thus, 
deep in the cool cores, the estimate above gives kiloparsec- 
scale temperature fluctuations, rising to tens and even hun- 
dreds of kpc at larger distances from the centre. 

Interestingly, temperature fluctuations on 1 to 10 
kpc scales have been detected in coo l -core clusters 
llSimionescu et al.l I2001I : iFabian et all I2OO6I : ISanders et all 
l2010al l while in the bulk of the cluster gas and in non-cool- 
core (radio-halo) clusters, the scales app e ar to be larger, 
around 100 kpc (jMarkevitch et all l2003l : iMiUion fc AUenI 
I2OO9I ). Thus, we again find the observed physical conditions 
intriguingly close to the marginal stability conditions set by 



As we already pointed out in section [5. II various macroscopic 
instabilities that play an important part in plasma dynamics, in- 
cluding those due to plasma effects such as anisotropic viscosity 
and thermal conductivity (MVI, MTI, HBI) act on time scales 
roughly comparable with the turbulence and are slow compared 
to the microinstabilities: e.g., HBI in cluster cores is estim ated to 
have growth times of order 10^ Myr l lParrish et al.ll2009h . 
^ CoUisional scales are intermediate between turbulence and 
plasma microphysics: the collision times are ~ 0.04 Myr, 

~ 0.001 Myr, u~} 1 Myr (the latter is the typical time 
for Ti and Te to equalize); the mean free path is X^fp ~ 0.01 kpc, 
where we have taken reference values of rie = 0.01 cm""^ and T = 



1 keV, collision frequencies are ex nT "^/^ and A^^fp 



1 J.2 
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plasma microphysics. Nevertheless, we would like to con- 
clude on a cautious note: whether the plasma contrives to 
satisfy the lower bound (|41|l by smoothing the temperature 
gradients or by aligning them carefully across the magnetic 
field remains unclear and underscores the need for a de- 
tailed theory of the nonlinear saturation of the GTI and 
other plasma microinstabilities. Observationally, it would be 
fascinating to see if any evidence can be obtained of correla- 
tions between the magnetic field direction and temperature 
fluctuations — presumably not an impossible task if one 
combines radio observations of polariz ed synchrotron em is- 
sion and X-ray temperature maps fcf. iTavlor et al.ir2006l '). 
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APPENDIX A: PLASMA PRESSURE TENSOR 

We start with the general kinetic equation for the distribu- 
tion function of a plasma species: 

where e is the particle charge, c the speed of light, a = 
{e/m){E-\-uxB / c)—Au/dt and C[f] is the collision integral. 
The flow velocity u appears in the kinetic equation because 
V is the peculiar velocity. Since {e/mc){v x B) ■ df/dv = 
—Cldf/di}, where Q is the Larmor frequency and i9 is the 
phase angle of the particle's gyration around the magnetic- 
field line, equation (|A1|I can be rewritten as follows 

We can now express the plasma pressure tensor P = 
m J d'^v vvf using the following identity 



peculiar velocity. We get, therefore. 



Pij = p±5ij — {p± —p\\)bibj 



M, 



ijkl 



AO. 



dt 



n Pm, I ~\~ ^In Pm k ^ ^ m — Ckl 



(A6) 



where Cki = mj d'^vvkViC[f] and we have introduced the 
heat flux tensor Qmki ~ nT- J d'^vv^VkVif . 

So far we have made no approximations. As promised 
in section [2] we now calculate all terms in equation (|A6|I 
assuming that we can use a gyrotropic (independent of -&) 
distribution function. This amounts to setting up a pertur- 
bation theory in which to lowest order, equation ^A2\ gives 
a gyrotropic equilibrium distribution, Qd/o/di} — 0, and at 
the next order we have QdSf/di} ~ ■ ■ ■, where only /o ap- 
pears in the right-hand side. The assumptions we need to 
achieve such an expansion are uj/Q <^ 1 and kp <^ 1 for all 
quantities involved. 

Since /o is gyrotropic, we may gyroaverage (vkvi) = 
(l/27r) / di3vkVi = {vl/2){5ki - bkbi) + vpkbi inside aU the 
velocity integrals in the square brackets in equation (IA6|) . 
so we get 



0, 



(A7) 



MijklPklV mUm — MijklCkl 
MijklPml'^mUk ~ 

-p± [b X (Vu) ■ (I + 3bb) - (I + 3bb) ■ {Vuf X b] , (A8) 

MijhlPmkV mUl — 

~pi_ [b X (Vn)^ ■ (I + 3bb) - (I + 3bb) ■ (Vu) x bj 

-4(px - p||) [b {b-Vuxb) + {b-Vux b) b] , (A9) 
dPki 



M, 



ijkl ' 



dt 



-4(Pi 



db ^ db ^ ^ 
b— xb+— xbb 
dt dt 



(AlO) 



Similarly gyroaveraging {vmVkVi) in the heat flux integral, 
Qmkl = q± {bmSki + 5mkbi + Smibk) — (3g± — q\\)hmhkbi. 
Therefore, 



uVmQmki = (I + 3bb) ■ [Vq^ + (Vq^)^] x b 
bx [Vq^ (Vq^)^] ■ (l-f-3bb) 

4(3gx - gii )^b {b ■ Vb x b) + {b ■ Vb x b) b] . (All) 



Assembling equations (|A7HA10)) and (|A11|I together in equa- 
tion (|A6|I , we obtain equations ([SHI}. 



vv = ^(\-bb) + vpb+^, (A3) 

T = (v^ib+ (vx xb) + (vx X b) (v^b + 
or, in index notation, Tij — {1/ 'i)MijkiVkVi, where 
Mijki = {Sik + Sbibk) ejinbn + mnb„ (Sjk + Sbjbk) . (A4) 
Therefore, after integration by parts with respect to i9, 

Pij =p±Sij - {p± -p\\)bibj - Mliti I d^vmvkVi (A5) 



We now substitute equation (|A2|) into the above expression 
and notice that J d?vvva ■ df/dv = after integration by 
parts and using the fact that J d?v vf = by definition of 



APPENDIX B: CGL EQUATIONS 

In order to derive equations ((3} and Q, we average equa- 
tion ()A2p over the gyroangles, (l/27r) J di?, which eliminates 
the left-hand side. In the remainder, we assume that the 
lowest-order distribution function is gyrotropic and so can 
be written as f — f{t, r, v, wy ). The time and spatial deriva- 
tives in equation (|A2|I are taken at constant v, so, in order 
for the gyroaverage to commute with them, they have to be 
transformed into derivatives at constant v and W|| , a nontriv- 
ial step because un = i; ■ b{t,r): 



db f df 
+ d-t-Hd^, 



(Bl) 
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(V/)„ = (V/U„+(Vb)..(|^ 



(B2) 



APPENDIX C: FLUX FREEZING 



Using these formulae and also df/dv — {v/v)df /dv + 
bdf/dv\\ and {v) = wyb, (vv) = (ui/2)(l - bb) + uffbb, 
we find that the gyroaveraged equation (|A2[I is 



(bb : Vu 
(V- u 



^ ^ 2 



V dv 



a ■ b 

df 
dv« 



V dv dv\\ 



2 V d- 



= C[f]. 



(B3) 



Changing variables from to {v±,v^^) or where 

/I — v'^/2B, transforms this equation into forms that are 
perhaps more f amiliar from t he well-known Kinetic MHD 
approximation ()Kulsrudlll983[ l. 

Equations (|3]) and Q are obtained by taking the ^1/2 
and vjj moments of equation (|B3P and integrating by parts 
wherever opportune. The coUisional relaxation terms are 
easiest to calcul ate with a simp lified coUision al operator, 
e.g., Krook ()Snvder et al.l Il997f ) or Lorentz (|Rosin et al.l 
I2OI0I ). To complete the picture, it may be useful to men- 
tion here that in some cases, especially when the pressure 
anisotropy p± —p\\ is small compared to the pressures them- 
selves, it is convenient to replace equations Q and by 
equation (|33p determining the evolution of p± — and an 
equation for the total pressure p = {2/3)p± + (l/3)p|| or 
temperature T defined by p = nT. Using equations @ and 
(HI, we get 



3 dT 1 dn , 
2 at n at 



1 dB 2 1 dn\ „ , 
5d^-3nd^)-^"^'(^^) 



where q ~ qj^+ '7||/2. The first term is compressional heat- 
ing, the second viscous heating and the third the heat flux. 
While the same-species collisions do not affect the evolution 
of temperature (because of the energy and particle conserva- 
tion), we do have to add to the above equation a temperature 
equilibration term, —{Z/2)niVig{Ti — Tg) for ions and nega- 
tive of the same for electrons, where Vie is the ion-electron 
collision frequency (the ion-electron temperature equilibra- 
tion terms were omitted in equations ([HJ and (U) because 
the relaxation of the pressure anisotropy was the dominant 
coUisional effect there) . In situations where radiative cooling 
is important (as in the case of galaxy clusters discussed in 
section [5^ . the electron temperature equation should also 



have a cooling term 



A(Te), w here A is the cooling 



function (e.g., lTozzi fc Normanl|200l| 

Note that, in principle, since we kept the FLR terms 
in the pressure tensor, we should also have kept FLR cor- 
rections in the CGL equations. These arise from the FLR 
contribution to the heat flux — in the coUisional limit, it 
is the usual dia magn etic heat flux 5q = (dnv^i-^/AQ.) b x VT 
(see iBraginskiil [l965l ) . While the unperturbed part of these 
FLR terms is small compared to other macroscale terms, 
their perturbed part is comparable to the perturbed gyrovis- 
cous stress terms (fc • 5G in equation (O). In equation (|14p . 
the diamagnetic heat-flux terms are part of the perturbed 
pressures 5p± and Since the instabilities we study in 
this paper are Alfvenically polarized and so are indifferent 
to pressure perturbations, we do not need to calculate the 
diamagnetic heat fluxes and, therefore, omit them. 



The non-MHD terms in equation (I12p will still preserve the 
magnetic-fleld topology if the electric fleld can be expressed 
in the form E = — iteff x B/c + Vx, where x is an arbitrary 
scalar function and liefi is some effective velocity fleld into 
which the flux will be frozen. Consider equation (|lip . It is 
not hard to show that the electron pressure term is 



V ■ Pe = Vpx, - V ■ [bb{p±e - P\\e)] 
= Vx(p±e -P\\e) - {P±e ~P\\e) 

^Vp||e + (P±e -P\\e) -g- 



B 



bVb 



(CI) 



where V_l = (I — bb)- V and we have used V-b = —b-'VB/B. 
Since the flrst two terms are perpendicular to B, they can 
be represented as a vector product of some effective vector 
fleld with B. Therefore, introducing 



Ucff = Ue 



eripB^ 



Vx(p_Le - Pile) 



'{P±e - Pile) ( -I- b ■ Vb 



X B, 



(C2) 



where = n— (c/47rene) V x B, we flnd from equation 

Uoff X B Vp[|e P±e-p\\e'VB 

h, = —. [L,6) 

c erie erie B 

With this electric field, Faraday's law (|10|) becomes 

dB ^ . ^, cVne X Vp||e 

— = V X (ucff- X B) 5 ^ 

ot eni 



, , „P±e - Pile \ VB 

+ ( cV 5^ I X — , 



(C4) 



which is an equivalent form of equation (|12[) . Thus, the mag- 
netic fleld is frozen into the effective velocity fleld UcS ex- 
cept for two effects. The first of the two n on-flux-con s ervini 
terms in equation (\C4^ is the well-known iBiermannI (Il95^ 
battery (believed to be one of the mechanisms responsible 
for seeding the cosmic plasma with the initial magnetic fluc- 
tuations, subsequen tly amplified by turbulent dynamo; see 
iKulsrud et al.|[l997l ): the second is an effect due to the elec- 
tron pressure anisotropy. While it is not a battery in the 
sense of producing magnetic field from a zero initial condi- 
tion (plasma is unmagnetized when B = 0, so there is no 
pressure anisotropy), it is independent of the fleld strength 
and, therefore, can act as a source term. 

The flux unfreezing effect due to the electron pres- 
sure anisotropy is small except at very small scales. 
If we use equation (|33[l to estimate p±e — P\\e ~ 
(rieTe / / B)dB / dt, we flnd, very roughly, that the flux 
conservation is significantly violated only if the scale of vari- 
ation of Te and B perpendicular to B is < (peAmfp)^''^- 
While this is larger than the electron inertial or Lar- 
mor scales, where flux normally unfreezes in a coUisionless 
plasma, and can be larger than the Ohmic resistive scale, it 
is still extremely small compared to any scales relevant for 
the macroscopic dynamics (for the reference cluster core pa- 
rameters used in section [5?2] we get l±_ ~ lO'' km). Note that 
the parallel flrehose and gyrothermal instabilities considered 
in the main part of this paper are unaffected by this flux un- 
freezing effect because we set fc_L = and the instabilities 
contained no perturbation of the magnetic-fleld strength. 



